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We apply a recently developed renormalization group (RG) method to study synchronization 
in a one-dimensional chain of phase-coupled oscillators in the regime of weak randomness. The 
RG predicts how oscillators with randomly distributed frequencies and couplings form frequency- 
synchronized clusters. Although the RG was originally intended for strong randomness, i.e. for 
distributions with long tails, we find good agreement with numerical simulations even in the regime 
of weak randomness. We use the RG flow to derive how the correlation length scales with the 
width of the coupling distribution in the limit of large coupling. This leads to the identification of a 
universality class of distributions with the same critical exponent v. We also find universal scaling 
for small coupling. Finally, we show that the RG flow is characterized by a universal approach to 
the unsynchronized fixed point, which provides physical insight into low-frequency clusters. 
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I. INTRODUCTION 

There has been much interest in the spontaneous syn- 
chronization of nonlinear oscillators [l|, [2], |3|. Synchro- 
nization is the phenomenon of how oscillators with dif- 
ferent intrinsic frequencies can oscillate at the same fre- 
quency due to the interaction between them. As the 
coupling between them increases, the system becomes 
more synchronized, i.e. a given oscillator is more likely 
to be entrained with others. This phenomenon is found 
in many different contexts, including Josephson junctions 
[H , lasers Q , neural networks Jp] , chemical oscillators []J , 
and even rhythmic applause [7|. 

Although systems that exhibit synchronization usu- 
ally involve driven and dissipative dynamics, they may 
be analyzed using ideas from equilibrium statistical me- 
chanics. For instance, a system can experience a phase 
transition from the unsynchronized to the synchronized 
state, in which a macroscopic fraction of the oscillators 
follows the same frequency. This is commonly known 
as the entrainment transition [H, Q and is analogous to 
the paramagnetic-ferromagnetic transition in spin mod- 
els [To] - The coupling between oscillators acts like the 
inverse temperature: as coupling increases, the system 
becomes more ordered. 

In the mean-field case, when each oscillator is coupled 
to all others, the entrainment transition occurs at a finite 
value of the coupling [l], Q . One may also consider finite- 
dimensional lattices of oscillators, and it has been shown 
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that the lower critical dimension for entrainment is two 

0, EH) E3- For d < 2, macroscopic entrainment exists 
only at infinite coupling, but for finite coup ling , there 
are still local frequency-synchronized clusters [13j |. 

A previous work presented a real-space renormaliza- 
tion group (RG) approach for a one-dimensional chain 
of oscillators [14| . It was successful at predicting cluster 
structure and frequency in the case of strong randomness, 

1. e. when the frequency and coupling distributions had 
long tails. It was thought that strong randomness was 
required in order for the perturbative decimation steps 
to be accurate. The numerical RG was successful over 
a wide range of distribution widths and accurately pre- 
dicted the scaling of the correlation length with distri- 
bution widths. An advantage of this RG approach was 
that it allowed the couplings to be randomly distributed. 
While most studies have assumed a constant coupling 
across the lattice, a physically realistic system would have 
disorder in the coupling Q. 

In this paper, we extend the analogy between synchro- 
nization and equilibrium physics even further by demon- 
strating universal features in the one-dimensional model. 
Universality is the existence of a class of systems that 
exhibit the same scaling behavior. For example, in spin 
models, the scaling of the correlation length with tem- 
perature near a critical point may be identical across dif- 
ferent physical Hamiltonians. In general, the presence 
of universal behavior is important because it means that 
the physical properties being studied are in some ways 
independent of the microscopic details. In equilibrium 
physics, there is a close relationship between universal- 
ity and the RG: the universal properties near a phase 
transition are explained by the RG flow there. 

Here, we find a similar relationship between universal- 
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ity and the RG in the Id nonequilibrium synchronization 
problem. We first show that the RG is accurate over a 
wide range of frequency and coupling distributions, even 
in the regime of weak randomness. Then by studying 
the flow of the RG for large and small coupling, we find 
how the correlation length scales with the width of the 
coupling distribution. This leads to the identification of 
several distinct universality classes, based on generic fea- 
tures of the frequency and coupling distributions. In fact, 
universality exists even far from the synchronized unsta- 
ble fixed point. By studying the flow towards the unsyn- 
chronized stable fixed point, we find a universal approach 
to it, which dictates the dynamics of low- frequency clus- 
ters independently of initial disorder realizations. 

This paper is organized as follows. Section II reviews 
the RG and reports on its performance for weak random- 
ness. In Sees. Ill and IV, we study the correlation length 
for large and small coupling, respectively. Section V ex- 
amines the universal approach of the RG to the stable 
fixed point. In Sec. VI, we conclude. 

II. RENORMALIZATION GROUP 

A. Overview 

The one-dimensional chain of oscillators with nearest- 
neighbor interactions is described by the equations of mo- 
tion: 

h = uj t + Ki-i sin(0i_i - 9 t ) + K t sin(0 i+ i - 9 t ) , (1) 

where 9i is the phase of the ith oscillator. The u>i are the 
intrinsic frequencies taken from a random distribution 
, assumed to have zero mean without loss of generality. 
The Ki are the couplings drawn from a random distribu- 
tion pk and are assumed to be positive. The couplings 
organize the oscillators into clusters of common frequency 
u) 1 defined as 

hm *<(*)- W. (2) 
(t-io)— >oo t — to 

In the renormalization group [141 ] . oscillators are com- 
bined into effective oscillators, so the model is slightly 
generalized to: 

rriiOi =rriiWi+Ki-ism(6i-i—9i)+Kism.(6i + i—9i) , (3) 

where the parameter rrii represents the number of original 
oscillators in the effective oscillator i. mi is referred to 
as an oscillator's mass, and the reduced mass between a 
pair of oscillators is given by a = m ', mi+1 . 

The RG is based on two decimation steps, which cor- 
respond to a coarse-graining of the system. The first is 
the strong coupling decimation step. Two oscillators, n 
and n + 1, connected by a large coupling K n would be 
expected to synchronize. Hence they are combined into 
a single effective oscillator with mass M, 

M = m n + m„+i , (4) 



and intrinsic frequency f2, 

O = {m n uj n + m n+1 uj n+ i)/M . (5) 

This step is valid up to zeroth order in ratios such as 
LOn/Kn and K n -i/K n . We refer to such a decimated 
pair of oscillators as strongly coupled. The new oscilla- 
tor may continue to be recombined with other oscillators, 
and hence oscillators that are strongly coupled belong to 
the same frequency-synchronized cluster. For the sake 
of clarity, when a pair of oscillators are decimated as 
strongly coupled, we say that their bond has been deci- 
mated. 

The second decimation step is the fast oscillator deci- 
mation step. An oscillator with large intrinsic frequency 
to n relative to its neighbors is expected to rotate freely. 
Such an oscillator will not synchronize with its neigh- 
bors, so the coupling to the neighbors, K n —i and K ni 
are set to zero. This step is valid up to first order in 
ratios such as K n /uj n and K n -\/u) n . In practice, the fast 
oscillator is removed from the chain and stored for later 
analysis, and the coupling between neighbors n — 1 and 
n + 1 is set to zero. There is a second-order shift in in- 
trinsic frequency for the fast oscillator and its neighbors 
[3], but we ignore the shift in this paper because it is 
small. The fast oscillator may consist of multiple oscilla- 
tors that were strongly coupled together. Hence, removal 
of a fast oscillator means that the cluster is in its final 
form. 

The chain of oscillators is renormalized by succes- 
sive application of the two decimation steps. The pro- 
cess is carried out numerically on a list of parame- 
ters {mi,u>i, Ki}. To decide which oscillator or bond 
to decimate first, the energies of each are calculated: 
{\u>i\, 2 t _! <i +1 }• Energy in this context reflects how much 
influence something has on clustering. An oscillator with 
the highest energy in the chain has a large frequency and 
should probably be decimated as a fast oscillator. Sim- 
ilarly, a bond with the highest energy should probably 
be strongly coupled. The highest energies are decimated 
first, so that coarse-graining corresponds to decreasing 
the energy scale of the system. In practice, we decrement 
the energy scale E and consider oscillators and bonds 
with energies > E for decimation on each cycle. 

Although an oscillator or bond may have the high- 
est energy, it must satisfy another criterion involving its 
neighbors before being decimated. We calculate the ratio 
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Pn,n+l\U n — LO n +l\ 

which measures the tendency of a pair of oscillators to 
synchronize. If Ip 1 > E and r n > 1, the bond is strongly 
coupled. If \u n \ > E, r„_i < 1, and r n < 1, oscillator 
n is removed as a fast oscillator. The threshold of 1 was 
chosen based on a numerical study of small chains [l5[ 
and is exact for a chain of two oscillators. 

It is possible that an oscillator or bond is not immedi- 
ately decimated despite having the largest energy in the 
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chain. For instance, if r n _i < 1 but r n > 1, oscillator 
n is not decimated as fast. Eventually though, the en- 
tire chain is decimated, leaving a list of fast oscillators. 
Each of these fast oscillators corresponds to a frequency- 
synchronized cluster. One may then study the statistics 
of cluster mass and frequency. 

The clusters predicted by the RG are compared to 
those found by numerically integrating Eqs. {l} with a 
variable stepping Runge-Kutta algorithm [16!] . To iden- 
tify frequency clusters in the simulations, the average 
frequencies are calculated according to Eq. A group 
of oscillators is determined to be a synchronized cluster 
if its members have the same value of uj within some 
tolerance. 



TABLE I: The critical exponent v, calculated by simulation, 
numerical RG, and analytical RG for different distribution 
types for ui and K. The first two columns list the distribution 
types: rectangular, triangular, Gaussian, Lorentzian, or delta 
function. The exponent v describes how the correlation length 
scales with the coupling width: £ ~ A". 

Puj Pk Simulation Numerical RG Analytical RG 



rec 


rec 


0.671(8) 


0.666(3) 


2/3 


tri 


tri 


0.673(8) 


0.668(2) 


2/3 


gau 


gau 


0.69(1) 


0.669(2) 


2/3 


lor 


lor 


0.47(1) 


0.48(2) 


1/2 


lor 


rec 


0.51(2) 


0.503(3) 


1/2 


rec 


lor 


0.668(8) 


0.669(3) 


2/3 


gau 
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2.04(5) 


2.003(4) 
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B. Application to weak randomness 

In the previous work [lj], the RG was found to be in 
good agreement with simulations in the regime of strong 
randomness. Comparisons of real-space cluster structure, 
cluster mass distribution, and cluster frequency distribu- 
tion were excellent. Distributions of both frequency (p u ) 
and coupling (px) were assumed to be Lorentzian. The 
long tails of the Lorentzian distributions heightened the 
overall randomness in the chain, helping to ensure accu- 
racy of the perturbative decimation steps and to support 
the intuitive notion of the strong randomness case. For 
instance, when decimating oscillator n as fast, uj n should 
be much larger than the neighboring lo, K . 

To understand the validity range of the RG, we re- 
lax the condition of strong randomness by considering 
cases where p^, and pk have finite variance. In particu- 
lar, we use rectangular, triangular, and Gaussian distri- 
butions. We still find good agreement with simulations in 
the regime of weak randomness. Figure[T]gives real-space 
comparisons between RG and simulation for rectangular 
and triangular distributions. Figure [2] compares RG and 
simulation in terms of cluster mass distribution. This 
shows that the RG is applicable over a wide range of 
distribution types and widths. 

The numerics were done with chains of 10 6 and 10 5 
oscillators for RG and simulation, respectively, p^ is as- 
sumed to be symmetric around zero, while pk is the pos- 
itive half of an otherwise symmetric distribution. The 
width of p u is defined to be the half width at half max- 
imum for Lorentzian, triangular, and rectangular, and 
the standard deviation for Gaussian. The width of px 
is defined similarly. The width of p^ is set to one, with- 
out loss of generality, while the width A of pk is varied. 
The RG is in good agreement with simulations for a wide 
range of A. In Sec. |Vl we provide insight into why the 
RG works so well even for weak randomness. 

With the simulation data, one can study the shape 
of the mass distribution p m . For large coupling width 
(A > 2.5), it takes the form 

p m (m) ~ m Cl e~ m ^ , (7) 



where C\ and c-i are constants, and C\ is approximately 
2/3 whenever p^ and px are of weak randomness. For 
smaller values of A, the mass distribution has a similar 
but more complicated form. 

In the limit of large coupling width, the correlation 
length £ of the system may be defined as the average 
cluster mass. In Fig. [3j we plot £ as a function of A for 
rectangular and triangular distributions. For all distri- 
butions considered, the scaling 

£ ~ r (8) 

holds for large A, where v is the critical exponent. Table 
U lists the values of v for different types of distributions. 
In Sec. Ill C| we comment on the significance of v. Note 
that A is a proxy for coupling strength, since larger A 
implies larger couplings in the system. 

C. Fixed points 

The RG described above is a functional RG since the 
decimation steps modify distributions of uj, K, and m, as 
opposed to a small set of parameters in a uniform system 
|17| . The flow of the RG is given by how those distribu- 
tions evolve as energy decreases. This corresponds to a 
coarse graining due to the buildup and removal of clusters 
through the decimation steps. Since there is not global 
synchronization in one dimension [l3l ]. clusters cannot 
grow forever. Indeed, by the end of the RG, the chain 
has been completely divided into a list of fast oscillators. 

The fixed points of the functional RG are the distribu- 
tions of uj and K that are invariant under the flow of the 
RG. We identify two fixed points of the RG. The stable 
fixed point is given by A = 0, or pk(K) = S(K). This 
corresponds to the unsynchronizcd phase, since all the 
oscillators are freely rotating. The unstable fixed point 
is given by A — oo. This corresponds to the synchronized 
phase, since p u is then relatively narrowly peaked at zero. 
Any finite A will flow to the unsynchronized fixed point 
(Figs. 2] and [5]). See Sec. IV Al for further discussion of the 
unsynchronized fixed point. 
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FIG. 1: Average frequency along the chain when u and K are both drawn from (a) rectangular and (b) triangular distributions. 
The RG predictions (open circles) are compared with simulation results (dashed lines, solid squares). The coupling width A is 
(a) 1 and (b) 10. 
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FIG. 2: Cluster mass distribution when lo and K are both drawn from (a) rectangular and (b) triangular distributions. The RG 
predictions (open circles) are compared with simulation results (solid squares). The plots for different values of the coupling 
width A are offset for visibility. 




FIG. 3: Correlation length £ vs. coupling width A when lo and K are both drawn from (a) rectangular and (b) triangular 
distributions. The RG predictions (open circles) are compared with simulation results (solid squares). The lines show the 
power law fits to the solid squares with the exponent v = 0.67. Here, £ is defined as the average cluster mass. 
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FIG. 4: Flow diagram of the RG, showing the unsynchronized 
stable fixed point at A = and synchronized unstable fixed 
point at A = co. 



The system may be interpreted as having a phase 
transition at A = oo, similar to the transition the one- 
dimensional Ising model has at zero temperature. The 
coupling width here is analogous to inverse temperature: 
the system becomes more ordered as A increases. The 
scaling in Eq. (JS} for large A probes the criticality of the 
system, and the correlation length £ diverges at the crit- 
ical point. 

It is common for equilibrium systems to exhibit uni- 
versality near the critical point. This means that macro- 
scopic quantities, such as the correlation length, scale 
with temperature in the same way for different systems 
near the critical temperature Although the model 

given by Eq. ([1} is far from equilibrium, since it is driven 
and overdamped, one may still find universal behavior 
near the critical point. Indeed, Table Q] shows that the 
exponent v is the same across different types of disorder. 
When and px are rectangular, triangular, or Gaus- 
sian, v « 2/3, and thus those systems belong to the same 
universality class. On the other hand, the Lorentzian 
case has v ss 1/2 and is in a different universality class. 

The significance of universality is twofold. First, it 
means that the scaling given by v does not depend on the 
microscopic details of a particular system. This is impor- 
tant when designing experiments to study synchroniza- 
tion. When it comes to observing the scaling in Eq. {SJ, 
the exact shapes of p u and px are not critical, as long as 
they are within a given universality class. Second, from a 
theoretical standpoint, universality is important in that 
it is often closely related to the RG flow near the critical 
point [l0(. In Sec. IIII1 we draw the connection between 
universality and the RG flow. The fact that the RG cor- 
rectly predicts the scaling laws boosts the claim that it 
is in fact a good representation of the model. 



the coupling width A in the regime of large A. The expo- 
nent v is seen to depend on generic features of the distri- 
butions of u> and K. Systems with the same such features 
will be in the same universality class. The approach here 
is similar to the usual one in equilibrium statisical me- 
chanics, i.e. examining the RG flow near criticality. 
A. RG flow 



As explained in Sec. Ill Al oscillators and bonds are 
decimated in order of decreasing energy. To simplify the 
analysis in this section, we define the energy of a bond to 
be K instead of K/2p. This means that a bond may be 
considered for strong coupling earlier on in the RG. This 
makes a difference only when the bond's r = K/h\Aoj\ > 
1 when E = K but would have been r < 1 when E = 
K/2p. In practice, this happens rarely. Empirically, it is 
more accurate to use K/2pL as the energy, but using K 
instead overestimates £ by only up to 5%. 

When A 3> 1, most bonds have higher energy than all 
the oscillators (Fig. [5]). Also, most bonds will satisfy r > 
1 and will be strongly coupled. Hence, the initial stage 
of the RG (when E > E* for some E*) is given by the 
strong coupling of all bonds that satisfy K > E. Since 
the strong coupling decimation is done independently of 
the m and ui of the corresponding oscillators, there are no 
correlations between oscillators remaining in the chain. 
E* is defined as the energy scale at which the presence 
of fast oscillators becomes important (when r ~ 1). In 
the remainder of this section, we calculate the RG flow 
for E > E*. 

Let p m (m, E) and px(K,E) be the normalized dis- 
tributions of m and K of the chain at a given en- 
ergy E. As the energy is decremented from E to 
E — dE, dE ■ pk{E, E) pairs of oscillators are decimated 
as strongly coupled, and the corresponding couplings dis- 
appear. The flow of pk(K,E) is therefore given by sim- 
ply rescaling on the interval < K < E: 



IIS w\ Pk(K,E ) 
Pk(K,E) = 

J h dK'p K (K',E ) 



(9) 



III. UNIVERSALITY FOR LARGE COUPLING 

In this section, we provide an analytical understanding 
of the power-law scaling of the correlation length £ with 

I 



where Eq is the initial energy and pk{K, Eq) is the initial 
K distribution. The flow of p m is given by successively 
combining pairs of oscillators: 



p m {m, E - dE) 



p m {m,E) + dE p K (E,E) [-2p m (m,E) + J J dm 1 dm 2 p m {m ll E)p m (m 2l E)6(m - (mi + m 2 ))] 

l-dEp K {E,E) 
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which leads to the integro-differential equation: 

dp. 



dE 



Pk{E,E) 



p m (m,E) 



dm'pm(m', E)p m (m — m',E) 



We treat to as a continuous variable on the range [ 0, oo). 
The initial condition at Eq can be approximated by 

p m (m,E ) = e- m , (12) 

which captures the fact that clusters have size one. The 
equation may be solved by Laplace transforming with 
respect to to. The solution is 

p m(m ,£)= * e-W^), (13) 

where 

t{E) = e -lE dE'pK(E' : E') ^ 



J E dK'p K (K',E ) 

is the average cluster mass at energy E. According to 
Eq. (|13[) . to is exponentially distributed in the initial 
stage of E > E*. (The final distribution of cluster 
mass is not strictly exponential, as seen in Eq. (J7JI and 
Fig. [2j) Equation (fT5|) relates energy to length scale. If 
Pk{K 1 Eq) is continuous and nonzero near K = 0, 

i ~ 4 when E < A. (16) 
E 

This is the case for triangular, Gaussian, rectangular, 
and Lorentzian distributions, although the proportional- 
ity constants differ. 



The £ and ao-dependence can be removed from the inte- 
grand by scaling the integration variables appropriately, 
and the resulting integral can be done numerically: 

Ou|Aw|) = 0.42 (7 £ 1/2 . (20) 

The length scale at E = E* is given by £(E*) = £. Com- 
bining Eqs. dHJ), (QUI, and (H, we find 

E* ~ e /2 ~ \±) ■ (21) 

Solving Eq. flU]) for E* yields 

- A 1/3 (22) 



(11) 

I 

B. Correlation length 



To calculate the correlation length £, we find the energy 
scale E* at which the system looks unsynchronized. It is 
given by the point in the RG when (r) ~ 1, since that is 
when the RG starts to encounter a lot of fast oscillators, 
which end cluster formation. Hence, 



E* ~ ( M |Au,|> , (17) 



where the expectation value is over all pairs of oscillators 
at energy E* . Since there are no correlations between 
to and u of different oscillators, we can use Eq. (fT3")) for 
both oscillators in the pair. Also, an oscillator of mass to 
has an to given by the average of to original frequencies. 
When to 3> 1 , we can apply the central limit theorem, to 
find the distribution 



p m .Am^E) = Q e -/«)(^r*#) (18) 



where ao is the standard deviation of the original p t 
Thus the expectation value is given by 



(19) 

I 

and thereby 

f - A 2/3 , (23) 

which predicts that v = 2/3 and matches simulation re- 
sults (Table H]). Although the exponent is the important 
quantity, we note that when all proportionality constants 
are included, the value of £ predicted here is within 10% 
of the simulated result. We also note that E* is universal, 
according to Eq. ([22]) . 

The universality of v among Gaussian, triangular, and 
rectangular distributions is seen to come from Eqs. (|16|) 
and (|18p . In other words, if p u has finite variance and 



( M |A„|) = 



dmidm 2 dojidoj2 p m , u (mi, uji, E) p m ,u(m2, w 2 , E) 



TO4TO2 

TOl + TO 2 



\UJ\ - L0 2 \ 
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C. Application to other distributions 



to 



Pk 

■ i 

i 



1 



X E 



FIG. 5: An example of when the initial coupling distribution 
(dashed-dotted line) is wider than the initial frequency distri- 
bution (dashed line), i.e. when the coupling width A > 1. 
Rectangular distributions are shown here. The RG does 
mostly strong coupling decimation until some energy E* , 
which corresponds to a length scale £ at which the system 
looks unsynchronized due to the emergence of fast oscillators. 



px is continuous and nonzero at K = 0, then v — 2/3. 

Now we check when the above argument is self- 
consistent. We assumed that A is large, so that the RG 
does only strong coupling decimation in the beginning. 
Also, £ should be large so that it is valid to treat m as 
a continuous variable, to use the central limit theorem, 
and so that Eq. (fTo)) would hold at E* . To find how large 
A should be for the power-law scaling to be accurate, we 
require £ ^> 1 in Eq. (|23]) in a self-consistent way. For 
rectangular p^ and px, we find the condition A 3> 0.24. 
According to simulation results in Fig. [3J the power-law 
scaling holds for A > 7.5. 

One may use the above results to calculate the dy- 
namical exponent z for the system. It describes how the 
diverging time scale of the system scales with the diverg- 
ing correlation length near criticality [l8|. Intuitively, 
when the chain is composed of long clusters, the clusters 
tend to have small frequencies, so the characteristic time 
scale of the system is large. The characteristic frequency 
of the final list of clusters may be approximated by the 
frequency standard deviation of the original oscillators at 
E*: 



J J dmduj p 7rltUJ (m,uj, E*) mu> 2 
J dm p m (m, E*) m 



(24) 



The expectation value is done over the original oscilla- 
tors, as opposed to clusters of oscillators, which is why 
the integrands include the factor m. In the case when 
Eq. (fT8"|) holds, the time scale is given by 



1 



1/2 



(25) 



Thus, the dynamical exponent z = 1/2 when p u has finite 
variance. This matches well with simulation results. 



When p u is Lorcntzian, a® is infinite, so we cannot 
use the central limit theorem to characterize oscillator 
frequencies as in Eq. (fT5|) . Instead, we use the fact that 
the average of random variables drawn from a Lorentzian 
distribution is described by the same Lorentzian distri- 
bution. Thus, Eq. (|18j) is modified to 



p m<u (m,u),E) 



1 

r 



/ e )p u (u,Eo) , (26) 



so that the second factor is independent of m. Fol- 
lowing the same procedure as before leads to the expo- 
nent v = 1/2, which agrees well with simulations (Table 
[3J. Note that v = 2/3 if p u is rectangular and pk is 
Lorentzian. 

We can extend the analysis to the case where all the 
couplings are equal to a constant, here denoted by A. 
This is the most commonly studied case in the literature. 
First we consider the situation when px is a rectangle of 
width e at K — A: 



Pk(K,E ) 
Equation ([To]) gives 
i{E) = 



A - e < K < A 



A-e < E < A 



E- A + e 

If Pu has finite variance, then Eq. (f2"Tj) becomes 

1/2 



E* 



£l/2 = 



E* - A 



(27) 



(28) 



(29) 



which, in the limit of e — > 0, leads to v = 2, which is 
close to the simulated value 2.04(5). This agrees with 
the value given in Q , based on the linear approximation 
to Eq. ((T|). The critical exponent can also be predicted 
using the result that for a chain of size N, the critical 
cou plin g for complete synchronization scales as 0(y/N) 
[HHi]. This implies that for an infinite chain, the length 
scale of synchronized clusters scales as 0(A 2 ). 

To summarize, we have derived the scaling behavior 
near the synchronized fixed point by studying the RG 
flow. The analytical predictions for the critical exponent 
match closely with the results of simulations and numeri- 
cal RG (Tablc|T|. We have presented a general procedure 
for handling different types of distributions, and it may 
be applied beyond the cases considered here. For exam- 
ple, one may study the case where px diverges at K = 0, 
which would require the scaling in Eq. (|16p to be modi- 
fied. 



IV. UNIVERSALITY FOR SMALL COUPLING 

In this section, we calculate the correlation length £ 
for small coupling and show that universality also exists 
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in this regime. This is surprising, since universality is 
usually found in the vicinity of a critical point. This 
discussion generalizes the estimate for £ given in Ref. [3 
to non-Lorentzian distributions and is cast from an RG 
point of view. 

Since the average cluster mass approaches one in the 
limit of A — > 0, we use a different definition for the corre- 
lation length based on the exponential decay of the final 
mass distribution p m (m): 



Pm(2) 
Pm(l) 



(30) 



It is not clear if p m has the same exponential form for all 
m, because it is difficult to collect simulation statistics 
for clusters made of more than two oscillators for A <C 1. 
In this regime, there are very few such clusters. Hence, 
the above definition is restricted to m = 1, 2. Now, £ can 
be less than one and in fact approaches zero when A — > 0. 

The RG flow can be used to calculate £(A) as in the 
previous section. Since px is now much narrower than 
p u at the beginning of the RG, most of the oscillators will 
be decimated as fast oscillators. Only when E < A will 
bonds start to be strongly coupled. If a bond satisfies 
r = K/n\Au>\ > 1 initially, it will still be so when the 
energy is at the bond's energy, because removing neigh- 
boring fast oscillators does not affect this bond's r. Since 
a cluster of mass m requires a consecutive sequence of 
m — 1 strongly coupled bonds, e -1 ^ in Eq. (|3"0|) is equal 
to the probability of r > 1 given by the initial frequency 
and coupling distributions: 



-i/C 



2 A' 



dK PK {K) / dA pa(A) , (31) 



where we have set /i = 1/2, and 



Pa(A) 



du p u (u) p u (u - A) A>0 (32) 



is the distribution of the absolute value of frequency dif- 
ferences between neighboring oscillators. 

If pk does not have long tails, the integral over A is 
approximately 2Kp/±(Q). This gives 



-Vf = A(K) I dtupjcu) 
— cA , 



(33) 
(34) 



where the proportionality constant c depends on the dis- 
tribution types. The correlation length is then: 



l 



log A + log c 
1 



log A 



(35) 
(36) 



For very small A, the constant c drops out, and the form 
of £(A) is universal across Gaussian, triangular, and rect- 
angular distributions. 




FIG. 6: Correlation length £ vs. coupling width A for small 
A for Gaussian (circles), triangular (triangles), rectangular 
(squares), and Lorentzian (plus signs) distributions, calcu- 
lated from simulations. The prediction from analytical RG is 
plotted for weak randomness (solid line) and Lorentzian (dot- 
ted line). Lorentzian is in a different universality class from 
the others. Here, £ is defined according to Eq. (|30p . 



On the other hand, if Pk has long tails such that (K) 
is infinite, the approximation in Eq. (|33p does not hold, 
and Eq. (|31[) must be integrated directly. For example, 
when p^ and pk are both Lorentzian [14| . 



4A 



arctan K 



K 2 



(37) 



and the integral may be done numerically. 

Figure [6] shows that the predictions in Eqs. (|36| and 
([57)) agree well with simulations for the various distribu- 
tion types. It is clear that Lorentzian is in a different 
universality class for small A. Note that the universal- 
ity class is determined by whether (K) is finite: if p^ is 
Lorentzian while px is Gaussian, the scaling in Eq. (|36|) 
still holds. 



UNIVERSAL APPROACH TO THE STABLE 
FIXED POINT 



In this section, we examine the flow of the RG near the 
unsynchronized fixed point. The discussion here is differ- 
ent from that given in Sec. IIV1 where we assumed that 
the system started out already near the unsynchronized 
fixed point, i.e. coupling width A < 1. Here, we allow 
the system to start anywhere, including A 3> 1, and look 
at the flow of the RG after enough renormalization steps 
have been carried out so that the system approaches the 
unsynchronized fixed point. 

We return to using K/2p, = K as the bond energy 
instead of just K. 
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b) O O-O-O O (HK) O-O 

c) O O O-O O O O-O O O 

FIG. 7: As the RG progresses, the chain is punctured by 
more and more K = bonds. At the start of the RG (a), 
all oscillators are coupled to their neighbors. Near the end of 
the RG (c) , it is unlikely to have two nonzero bonds in a row. 
We study the properties of pairs of oscillators that are still 
connected with a nonzero bond. 



A. Unsynchronized stable fixed point 

At the unsynchronized fixed point, px{K) — 5(K), 
meaning that all the remaining oscillators are effectively 
uncoupled to each other. Since this fixed point is sta- 
ble, at late stages of the RG, the chain will be punc- 
tured by many K = bonds (Fig. [7]). This is due to 
the removal of a lot of fast oscillators, which leaves the 
former neighbors uncoupled. Due to the large fraction 
of K — bonds, most oscillators will be uncoupled to 
both neighbors. Such isolated oscillators are waiting to 
be removed as fast oscillators and will not be strongly 
coupled with their neighbors. Unless the u) distribution 
of isolated oscillators diverges at zero, it will approach 
a uniform distribution on the interval [-E, E] and zero 
elsewhere, since the flow is given by successively chopping 
off the sides. Thus at the unsynchronized fixed point, 
p u {u)) — j£ for \uj\ < E. Note that these isolated, low- 
frequency oscillators can be clusters of mass one or of 
higher mass, depending on the size of A. 

Due to the stability of the unsynchronized fixed point, 
almost all chains will flow to it, regardless of the initial 
distributions of u and K. Figures [5] shows examples of 
chains with different initial distributions flowing to the 
same fixed point. 

Now we consider the way p^ and p u approach their 
respective fixed distributions. We look at the component 
of the distributions that are not the fixed distributions. 
For pg, this means considering only bonds with K =/= 0. 
For p u , this means considering only oscillators with at 
least one nonzero bond. These components are pertur- 
bations that decay as the RG approaches the fixed point. 
Figures [HI show these perturbations. 

It is seen that the K perturbation collapses onto a 
triangular distribution regardless of the initial distribu- 
tions, while the u> perturbation collapses onto a flat and 
quadratic shape. Hence, there is a universal approach to 
the stable fixed point. One can think of the fixed point 
as having a least negative eigenvalue, corresponding to a 
favored direction for approaching the fixed point in the 
functional RG space. 



This suggests why the RG works well even with weak 
randomness. Regardless of whether the chain starts with 
strong or weak randomness, distributions of oscillator 
properties will look similar as the RG progresses. 

In the following section, we explain how the triangular 
K perturbation comes about. Due to rather technical 
details, we postpone the explanation of the flat-quadratic 
u> perturbation to the appendix. 



B. Explanation of the triangular perturbation 

Here, we explain how the RG algorithm leads to the 
universal shape of the K perturbation. Since the per- 
turbation is due to nonzero bonds, we consider pairs of 
oscillators with a nonzero bond connecting them, but 
otherwise uncoupled to the rest of the chain (Fig. [7J. 
We ignore single oscillators that are uncoupled to both 
neighbors, since they are waiting to be removed as fast 
oscillators and do not contribute to the buildup of future 
clusters. 

To facilitate the discussion, we define a center to be a 
pair of oscillators that both satisfy |w| < E, where E is 
the energy of the RG at a given stage. A tail is a pair 
of oscillators where at least one \u)\ > E. In other words, 
neither of the oscillators in a center has been checked for 
fast oscillator decimation. On the other hand, at least 
one of the oscillators in a tail has been checked for fast 
oscillator decimation, but since it is still in the chain, we 
know that it satisfies r > 1 and is waiting to be strongly 
coupled (when E — K). 

A pair can stop being a center by either being checked 
for fast oscillator decimation or strong coupling decima- 
tion. In the first case, the pair will either be removed 
(r < 1) or become a tail (r > 1). In the second case, 
the pair will either be strongly coupled (r > 1) or remain 
intact in the chain (r < 1). It can be shown that when 
a center gets checked for strong coupling decimation, it 
will always be decimated [201 ] . It is also possible for new 
centers to be formed, due to the strong coupling of other 
oscillators; however, we can neglect this possibility when 
it is unlikely to have two nonzero bonds in a row. 

Let the centers be described by the distribution 
p c (u>i, u>2, K; E), where |a>i| > 1 1 • This distribution is 
nonzero only when \u>i\, K < E. Its flow is given by set- 
ting to zero the region that doesn't satify that inequality, 
and then normalizing to unity. Unless there is a diver- 
gence at the origin, the distribution of centers approaches 

p c (u> 1 ,oj 2 ,K;E) = |wa| < |wi| < E,0 < K < E . 

This is just the statement that any smooth distribution 
will look flat if you keep chopping off its sides. Integrating 
out Ui, u>2 in Eq. ([55]) . we find that the distribution of K 
of centers is flat: 

pUK;E) = I 0<K<E. (39) 
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FIG. 8: Distribution of oscillator frequency (top) and bond energy (bottom) at different stages of the RG: (a,d) 50%, (b,e) 
25%, and (c,f) 3% of original oscillators left in chain. The results for different initial distributions of ui and K are shown: 
Lorentzian A = 7.5 (black, solid), triangular A = 1.25 (gray, solid), triangular A = 7.5 (black, dotted). The RG approaches the 
fixed point regardless of the initial chain: the distribution of uj becomes rectangular, while the distribution of K approaches a 
delta function at zero. 



We now consider centers at the moment they get 
checked for fast oscillator decimation, i.e. when uj\ = E. 
Here we assume that u>i is positive, but analogous results 
hold for when it is negative. 

pZ, i=E (u 2 ,K;E) = ^ 0<\u }2 \,K<E (40) 

At this point, if ui2 and K are such that r < 1, then 
oscillator 1 will get removed as a fast oscillator, and os- 
cillator 2 will become an isolated oscillator, which is then 
ignored. But if r > 1, then the pair becomes a tail and 
remains in the chain. Thus the new tails at energy E are 
described by: 

* , ~ s 1 1^2 - E\ 

pl i=E (u; 2 ,K:E) = — 0< 1 2 2 1 < K <E, 

(41) 

which can be rewritten as 

f& 1=E (\Au>\,K- t E) = ± 0< \M <K<E . 

. (42) 

This distribution is constant over a triangular region in 
\Au\, K as illustrated in Fig. [TO] 

Now here is the key. Tail pairs remain in the chain 
until they are decimated as strongly coupled, which hap- 
pens when the RG energy reaches the bond energy. Thus 



Eq. (|4"2"| describes the flow of the existing tails at any 
energy E. Existing tails have the same distribution as 
incoming tails at any energy, so Eq. (|4*2"|) describes all 
tails. 

By integrating out |Aw| in Eq. (|42|. one finds that the 
K distribution of tails is triangular: 



p\(K-E) = -^k 0<K<E. (43) 



In principle, the combined K distribution of center and 
tail pairs is given by a weighted average of Eqs. (|39| and 
(|43| . But it is shown in the appendix that as energy 
decreases, there are more and more tails relative to cen- 
ters, so that at low energy, all non-zero K are given by 
just Eq. (|43j) . This explains the triangular distribution 
in Fig. M 

The universality of Eq. (|43f comes from the applicabil- 
ity of Eq. (f3"5| . Unless diverges at ui = 0, will ap- 
proach 5(K) through a decaying triangular perturbation 
as energy decreases. This explains the universal approach 
to the stable fixed point. 
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0.5 1 1.5 0.5 1 1.5 0.5 1 1.5 

k/E k/E k/e 

FIG. 9: Distribution of frequency of unisolated oscillators (top) and nonzero bond energies (bottom) at different stages of the 
RG: (a,d) 50%, (b,e) 25%, and (c,f) 3% of original oscillators left in chain. The results for different initial distributions of uj 
and K are shown: Lorentzian A = 7.5 (black, solid), triangular A = 1.25 (gray, solid), triangular A = 7.5 (black, dotted). These 
are the perturbations that decay as the RG approaches the fixed point. The ui perturbation becomes quadratic in the middle 
and flat otherwise, while the K perturbation becomes triangular. This shows that the system approaches the fixed point in the 
same way, regardless of the initial chain. 
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FIG. 10: Illustration of the flow of tail pairs according to 
Eq. (|42[) . The distribution is constant on a triangular region 
in | Aw |, K space and zero elsewhere. As energy decreases, the 
triangular region shrinks, but keeps the same shape. Both ex- 
isting tails and incoming tails at a given energy are uniformly 
distributed on the same triangular region. 



C. Physical interpretation 

The above discussion has focused on RG flow at low 
energies, which corresponds to low frequency clusters. 



When E is sufficiently small, Eq. (|38[) accurately de- 
scribes centers. In physical terms, Eq. ([55)1 describes 
pairs of clusters with | u^i | , ] 0^2 ] < E that are interact- 
ing via a bond K. Depending on how big the coupling is, 
the two clusters may end up synchronizing and forming 
a bigger cluster. 

Most of the clusters with \uj\ < E will be effectively 
isolated from the rest of the chain, since most couplings 
are zero at low energy. It is shown in the appendix that 
of the cluster pairs that are coupled and have at least one 
cluster with \oj\ < E, 3/5 are centers and 2/5 are tails. 
The centers have interactions given by Eq. (j3T)|) , while 
tails have interactions given by (|43|) . Since Eqs. (|38|) and 
(|43[) and the proportion of centers and tails are universal 
across different initial distributions, we know that the 
dynamics at low frequency are also universal. Thus, this 
analysis of the RG provides insight into the interactions 
of low frequency clusters. 

A possible extension of this analysis is to study the dis- 
tribution of cluster mass. As mentioned in Sec. Ill Bl the 
final distribution of cluster mass is of the form in Eq. (J7J). 
In the current discussion, the mass information is stored 
in the bond energy K/2fi. We have considered the flow of 
pairs of oscillators that are uncoupled from the rest of the 
chain. By considering the flow of longer isolated chains, 
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one may be able to understand the buildup of larger clus- 
ters and hence the final mass distribution. The difficulty 
lies in the correlations that appear due to the r criterion. 

VI. CONCLUSION 

In this paper, we have explored various features of the 
real-space RG approach to Id synchronization, first pre- 
sented in Ref. [lj. We have shown that the RG method 
performs well even beyond the strong randomness case, 
for which it was originally intended. The RG was also 
used to calculate critical properties of random oscillator 
chains, such as correlation-length scaling for both large 
and small coupling. We identified several universality 
classes, whose behavior we also derived analytically. Ex- 
cellent agreement was found between our analytical ar- 
guments, numerical RG, and simulations. Finally, we 
demonstrated the universal approach to the stable fixed 
point. 

The results presented here may find relevance in phys- 
ical realizations of the one-dimensional model. The uni- 
versality implies that the predicted scaling could be ex- 
hibited in experimental realizations of the model with- 
out fine tuning: as long as the distributions of uj and 
K have the generic features of a given universality class, 
the corresponding scaling laws will hold. In addition, the 
dynamical exponent derived here shows the relationship 
between time and length scales in this model. 

We emphasize that the results were all based on the RG 
as opposed to a traditional dynamical-systems approach. 
A natural next step is to develop an RG for higher dimen- 
sional lattices. It has been determined numerically that 
the lower critical dimension for macroscopic entrainmcnt 
is two [1,0]. Thus for d > 2, an RG would have two sta- 
ble fixed points, corres pon ding to the synchronized and 
unsynchronized phases [l(| • Once an RG has been devel- 
oped for higher dimension, it may be possible to obtain 
an analytical understanding of synchronization in a way 
similar to the present work. By studying the RG flow 
near criticality, one may even find universal behavior in 
higher dimension. 

This work has been supported by Boeing. GR thanks 
the Research Corporation and the Packard foundation for 
their generous support. We also thank Heywood Tarn for 
many helpful discussions. 



APPENDIX A: FLOW OF TAILS AND CENTERS 

In this appendix, we calculate in more detail the flow 
of distributions of center and tail pairs. We consider 
late stages of the RG when it is unlikely to have two 
nonzero couplings in a row. We look at pairs of oscillators 
with a nonzero coupling between them, but otherwise 
uncoupled to the rest of the chain. Let lj% and u>2 be the 
frequencies of the pair, with |wi| > \uj2\- A center pair 
satisfies |wi|, \lo2\ < E, while a tail pair satisfies |wi| > E. 



Thus, at least one of the oscillators within a tail pair has 
been checked for fast oscillator decimation and failed. 
Since we are considering only center and tail pairs, once 
such a pair gets decimated as either strongly coupled or 
fast, it drops out of the discussion. 

We consider the rates of three processes as energy is 
decreased: the rate at which centers are decimated as 
strongly coupled, the conversion rate of centers into tails, 
and the rate at which tails are decimated as strongly cou- 
pled. In general, the distributions will not be normalized 
to one, since we would like to track the number density 
as energy is decreased. 

We repeat part of the discussion in Sec. IV Bl with gen- 
eralized notation for the sake of clarity. The starting 
energy is denoted by E, while the energy at a given stage 
of the RG is E. The energy at which a particular tail 
pair is created is denoted by Eq. Note that energies are 
always taken to be positive. 

1. Flow of centers 

We start the calculation by assuming that the distri- 
bution of all unisolated pairs is given by the fixed point 
distribution of centers, Eq. (|55|) . This means that all 
pairs are centers, and there are no tails. The number 
density of tails is 

P c {lo u uj2,K;E) = -t- \u 2 \ < |wi| < E,0 < K < E . 
2E A 

(Al) 

This distribution is normalized to (E/E) 3 . The number 
of centers decreases as E decreases due to the decimation 
of centers as strongly coupled or fast and the conversion 
of centers into tails. Note that we assume that no new 
centers are formed, which is a valid assumption during 
late stages of the RG when most couplings are zero. 

By integrating out variables in Eq. I|A1[) . one obtains 
the distributions of uj\ and u>2 ■ 

pZM^ E ) = il^il \u 1 \<E (A2) 

E- 3 

p%^ 2 ;E) = jL (E -\oj 2 \) \lj 2 \<E. (A3) 

Thus the number density of all center frequencies is 

fi(w;E) = pZ 1 (u;E)+pZ 2 (u;E) (A4) 

= £ H < E , (A5) 

which is normalized to 2(E/E) 3 , reflecting the fact that 
each center contributes two frequencies. 

2. Flow of tails 

Now we consider the transition from center to tail, i.e. 
when uji — E. We assume for now that uj\ is positive, 
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but analogous results hold for when it is negative. Im- 
mediately before checking the pair for fast oscillator dec- 
imation, the pair is described by 



p c (E Q ,m 2 ,K;E ) 



(A6) 



0<\lo 2 \,K <E a (A7) 



There is a half chance that the pair will satisfy r < 1 
and be decimated out, and there is a half chance that 
it will satisfy r > 1 and become a tail. Assuming the 
latter case, immediately after becoming a tail, the pair is 
described by 

E Q 2 

(A8) 

which reflects the fact that r > 1. As shown in Fig. [101 
the flow of tails created at Eq as energy E decreases is 
given by 



p\oj 2 ,K;E,E ) 



1 

El 



0< \^M <K <E 



(A9) 

By integrating out K above, one finds the distribution 
of oj 2 for tails that satisfy uj\ = Eq: 



pl+J^E^o) 



Eq — Ll>2 . 



(A10) 



En - 2E < Uo < En 



The corresponding distribution for tails that satisfy 
u>i = —Eq is 



pl 2 {uj 2 ;E,E ) = ^{E , 



(All) 



-E <lu 2 < -En + 2E . 



The distributions in Eqs. (|A10[) and (|A11[) are normal- 
ized to (E/En) 2 , so that the survival probability is one 
when E = Eq, but decreases as E decreases due to strong 
coupling decimation. Note that this probability can be 
deduced geometrically from Fig. [TO] 

We can also write down the distribution of uj\ for tails 
that were made at Eq: 



E,Eo) 



E 



(A12) 



— ) Sfa + Eo), (A13) 



for positive and negative respectively. 

Now we find /?* (u>;E), the distribution of all tail fre- 
quencies at a given energy E. We integrate over all the 
possible energies Eq, at which tails could have been cre- 
ated. As the creation energy is decremented from Eq to 
Eq — dEo, cIEq ■ p° (Eq; Eq) center pairs with positive uj\ 
are checked for fast oscillator decimation. An equal num- 
ber of center pairs with negative u>\ are checked. Note 
that p^ is given by Eq. (|A5[) . This leads to 



pU^e) 



E 



gIEq p c JEq;E ) [/£>; E, E ) + E, Eq) + p^w; E,E ) + ^(w; E, E ) ] (A14) 

plX (w; E) + j?~ (w; E) + pl+(uj; E) + ?~ (uj; E) . (A15) 



r 



The factor of 1/2 accounts for the fact that half of the 
pairs checked for fast oscillator decimation will survive 
and become tails. 

Then we calculate each of the terms: 



Pt + ME) 




rmin(E.u+2E) p 2 i 



\uj\ < E 
uj > E 



Eq 



(A18) 



Pl + M E ) = \ I dE §|j - Eq) 

2 J E E A Eq 

E < uj < E (A16) 
- E < uj < -E (A17) 



E 2 
2E 3 
E 2 
2E 3 



ffcME) 




Eq 



(A19) 



We have assumed that uj,E -C E. In other words, we are 
looking at low frequencies at low energies, i.e. when the 
RG is close to the unsynchronized fixed point. 
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Summing up the four terms, we find the frequency dis- 
tribution of tails in this regime: 

^'{JjfM^ ■ (A20> 

By comparing Eqs. (|A5|) and (|A20[) . one sees that there 
are relatively more and more tails than centers as E de- 
creases. By integrating p^(uj;E) and plj(uj;E) on the 
interval — E < uj < E, one finds that there are 2(E/E) 3 
and (2/3)(E/E) 3 oscillators with \oj\ < E belonging to 
tails and centers, respectively. In otherwords, there are 



(E/E) 3 center pairs and (2/3)(E/E) 3 tail pairs with at 
least one oscillator satisfying \ui\ < E. Thus, 3/5 of all 
such pairs are centers. 

By adding on the frequency distribution of centers 
given by Eq. (|A5[) . we find the frequency distribution 
of all unisolated oscillators 

KteB) - { JL H ■ (A21) 

This matches the flat-quadratic distribution found with 
the numerical RG in Fig. [9] 
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